Thermodynamical properties of dark energy 
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We have investigated the thermodynamical properties of dark energy. Assuming that the dark 
energy temperature T ^ a~" and considering that the volume of the Universe enveloped by the 
apparent horizon relates to the temperature, we have derived the dark energy entropy. For dark 
energy with constant equation of state w > — 1 and the generalized Chaplygin gas, the derived 
entropy can be positive and satisfy the entropy bound. The total entropy, including those of dark 
energy, the thermal radiation and the apparent horizon, satisfies the generalized second law of 
thermodynamics. However, for the phantom with constant equation of state, the positivity of 
entropy, the entropy bound, and the generalized second law cannot be satisfied simultaneously. 
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Results from numerous and complementary observa- 
tions show an emerging a paradigm 'concordance cos- 
mology' indicating that our universe is spatially flat and 
composed of about 70% dark energy (DE) and about 
25% dark matter. The weird DE is a major puzzle of 
physics now. Its nature and origin have been the intrigu- 
ing subject of discussions in the past years. The DE has 
been sought within a wide range of physical phenomena, 
including a cosmological constant, quintessence or an ex- 
otic field called phantom . Except the known fact that 
DE has a negative pressure causing the acceleration of 
the universe, its nature still remains a complete mystery. 
In the conceptual set up of the DE, one of the impor- 
tant questions concerns its thermodynamical properties. 
It is expected that the thermodynamical consideration 
might shed some light on the properties of DE and help 
us understand its nature. 

The topic on the DE entropy, temperature and their 
evolution by using the first law of thermodynamics was 
widely discussed in the literature 0, i, II H, H 0, B i, ■ 
It was found that the entropy of the phantom might be 
negative The existence of negative entropy of the 

phantom could be easily seen from the relation Ts — p+p 
between the temperature T, the entropy density s, the 
energy density p and the pressure p. Negative entropy is 
problematic if we accept that the entropy is in association 
with the measure of the number of microstates in statis- 
tical mechanics. The intuition of statistical mechanics 
requires that the entropy of all physical components to 
be positive. Besides if we consider the universe as a ther- 
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modynamical system, the total entropy of the universe 
including DE and dark matter should satisfy the sec- 
ond law of thermodynamics. The generalized second law 
(GSL) forphantom and non-phantom DE has been ex- 
plored in [8] . It was found that the GSL can be protected 
in the universe with DE. The GSL of the imiverse with 
DE has been investigated in [H, as well. In order to 
rescue the GSL of thermodynamics, Bekcnstcin conjec- 
tured that there exists an upper bound on the entropy 
for a weakly self-gravitating physical system . Beken- 
stein's entropy bound has received independent supports 
. A holographic entropy bound [l^ was subsequently 
built and it was argued to be a real conceptual change 
in our thinking about gravity jlj . The idea of the holo- 
graphic entropy bound was found to be a useful tool in 
studying cosmology (TBj . 

In the discussion of thermodynamical properties of the 
universe, it is usually assumed that the physical vol- 
ume and temperature of the universe are independent 
and by using the integrability condition d^S/dVdT = 
d^S/dTdV and the first law of thermodynamics, one ob- 
tains the constant co-moving entropy density. However, 
if we apply this treatment in the universe with DE, we 
found some problems of the DE thermodynamics 
Naively, we may think the DE temperature is equal or 
proportional to the horizon temperature Tu- It was 
found that the equation of state of the DE is uniq uely 
determined and the phantom entropy is negative [lOj. 
Therefore, a general DE model is not in thermal equi- 
librium with the Hawking radiation of the horizon. Be- 
sides, although the GSL can be valid for w > —1, for 
the phantom with w < —1, it was found that the GSL 
breaks down due to the negative temperature deduced 
in the formalism where the volume and the temperature 
are assumed to be independent [l0[. In summary, for 
the phantom, we either run into negative entropy prob- 
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lem or the GSL is violated. It is more realistic to con- 
sider that the physical volume and the temperature of 
the universe are related, since in the general situation 
they both depend on the scale factor a{t). In the cosmo- 
logical context, the apparent horizon is important, since 
on the apparent horizon there is the well-known corre- 
spondence between the first law of thermodynamics and 
the Einstein equation |16i] . On the other hand, it was 
found that the apparent horizon is a good boundary for 
keeping thermodynamical laws Q. Considering the ap- 
parent horizon as the physical boundary of the universe, 
it was found that both the temperature and entropy can 
be positive for the DE, including phantom. Furthermore, 
by considering the realistic case that the physical volume 
and the temperature are related, the GSL is proved to be 
always satisfied within the volume of the apparent hori- 
zon [lO|. Thus, in studying the DE thermodynamics, it 
is more appropriate to consider the universe in which the 
volume and the DE temperature are related. 

In this work we will investigate the thermodynamical 
properties of DE by assuming that the physical volume 
and the temperature are not independent. Now again it 
is natural to think that the DE is in thermal equilibrium 
with the Hawking radiation of the apparent horizon. In 
this case, we found that the DE entropy is the dominant 
entropy component and it becomes negative even for DE 
with w > —1 10]. Recall that the radiation temperature 



in the universe scales as T • 



, so we assume here that 



the DE temperature has a similar behavior T ^ a~" to 
avoid the negative entropy problem, where n is an ar- 
bitrary constant. It is not necessary to take n = 1 to 
ensure that the DE is in equilibrium with the thermal 
radiation, since their dispersion relations could be com- 
pletely different [1, From the above discussions, it 
is reasonable to expect that a physically acceptable en- 
tropy of the DE should be positive and satisfy the entropy 
bound. It should also satisfy the property required by the 
GSL. Since the usual thermal radiation temperature in 
the universe decreases as the universe expands, we expect 
that the DE temperature also preserves this property. 

By using the first law of thermodynamics TdS — 
dE + pdV for the DE, and considering the volume of 
the universe within the apparent horizon V — 47rf^/3, 
the total DE E — pV, we can express the DE entropy as 
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Pt^^^iPt + 3pt)dp, (1) 



where the Friedmann equation and the energy conserva- 
tion law have been used in the derivation, pt and pt de- 
note the total energy density and pressure, respectively. 
Taking derivative with respect to time on both sides of 
the above equation, we have 

S^27t(^^'^ '"\f'\p, + ip,)H{p + p)/T. (2) 

It can be seen that 5* > (< 0) if {p+p)/T > (< 0) dur- 
ing radiation dominated era (RD) and matter dominated 



era (MD), and S' < (> 0) if {p + p)/T>Q (< 0) during 
DE domination. In the phantom domination, the appar- 
ent horizon entropy decreases as the Universe expands 
[13]. This requires 5' > and ]! + 3w1Th/2 > T > to 
protect the GSL. Thus in the phantom domination era, 
the temperature of the phantom has to be positive to 
rescue the GSL. 

The radiation entropy can be obtained as usual Sr = 
sV , where s — a j a?" is the physical entropy density and a 
is the constant co-moving entropy. For DE with constant 
equation of state w, using the Friedmann equation, the 
entropies of the radiation and the apparent horizon arc 



Sr = Sr^X-'^^t(x)-'^l'^, 

Sa = Saq'^A^Y^ , 



(3) 



(4) 



where r2t(a;) = VLmQX^'^ -I- ilro^;"^ + ri^ox~'^(-'^+'^^ and 

To get the DE entropy we need to solve Eq. ([1]) by 
assuming T — Tol'j/ao) In the evolution of the uni- 
verse, the solution to Eq. ([T]) is given in the form 



Sw ^roTwQ 
SrQ flwoTrQ 



S: 
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(5) 

for the RD, MD and DE domination respectively, where 

Swl ^Ild Shj2 

are integration constants. 

As we mentioned previously, the intuition of the statis- 
tical mechanics requires positive entropy. We expect that 
this should also hold for the entropy of DE if it is sup- 
posed to keep the same microscopic meaning. From Eq. 
([5]) we learn that for DE with constant equation of state 
w > —1, non-negative Sw can be obtained if n > 3w — 3 
during RD. During DE domination, if n > —3(1 + w)/2, 
then Sw — oo when a — > oo. Thus, to get positive 
entropy for the DE, the parameter n should be chosen 
within the range —6 < 3w — 3 < n < —3(1 + w)/2 < 0. 

This parameter range of n can be further constrained 
if we express the solution of Eq. ^ as 



Sw 
SrO 



w 



^wQ TrO r 3+„_3 



^rO TwO 

a/ao 

(n — 3w) I X 



'n,{x)-^/' 
-^+^'^^n,{xr^/^dx 



(6) 



If n > 3w, the second term in the above equation is 
negative, which might lead Sw to be negative. There- 
fore, we need to restrict 3w — 3 < n < 3w to ensure the 
positivity of Sw- Note that for radiation, n = Sw = 1, 
Eq. ([6|) reduces to Eq. ([3]). Since n < 0, the dark 
energy temperature will increase with the scale factor 
a and at the present moment Two ^ TrO- During RD 
and MD, it can be seen that both DE entropy and the 
radiation entropy increase. However, if one notes that 
Sw/Sr = 3(1 -I- w)pwTr/i4:prTw) < 1, the DE entropy 
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FIG. 1: The evolution of for w = -0.9 and 5"^. The 
dotted hne is for {S^,/Sro) x {T^o x lO-^V^ro) with n = -5.0, 
the dash-dot line is for (Sw/Sro) x (r„o x 10~*/rro) with 
n — —3.5, the dash line is for Sr/Sro, and the solid line is for 
the apparent horizon entropy Sa/Sto- 

increases slower than the radiation entropy. During the 
DE domination, both the DE entropy and the radiation 
entropy decrease, 5*^, — > Sw2 > and Sr ^ when 
a oo, so > Sr in the future. Since the appar- 
ent horizon entropy increases during the DE domination, 
Sa = 3(1 -I- w)SaH, Syj = {n + 3(1 -t- w)/2)S^H and 
Sr — 3(1 + 3w)SrH/2, so the GSL is always respected 
for the DE with constant equation of state w > ~1. To 
see these points more clearly, we solve Eq. ^ numer- 
ically by choosing w = —0.9, flmo — 0.3, flwo — 0.7 
and firo — 8.35 x 10^^. The results for n — —5.0 and 
n = —3.5 are shown in Fig. [1] The numerical results 
confirm that by constraining iw — 3 < n < 3u', is 
positive. It is easy to see that the DE entropy and the 
radiation entropy are much smaller compared to the ap- 
parent horizon entropy, thus the entropy bound is always 
held. Although the radiation entropy and the DE entropy 
may decrease in the DE domination, due to their very 
small scale, their decreasing behaviors can be overcome 
by the increase of the entropy on the apparent horizon. 



Thus, including the total entropy in the universe and the 
entropy of the apparent horizon, we find that the GSL is 
protected. 

Now we come to consider the phantom with constant 
equation of state w < —\. In the RD and MD eras, 
if n < "iw — 3, is positive but it decreases starting 
from oo as the universe expands. The entropy bound 
is violated at very early times. Ifn > — 3(l-f?i;)/2, 
is negative during RD and MD and 5*^ cxd when 
a — > oo, so in the future although the GSL can be pro- 
tected, the entropy bound will be violated. If 3w — 3 < 
n < —3(1 + w)/2, Sw is negative. Thus, for the phantom 
with constant equation of state, it seems impossible to 
get a viable thermodynamics. The requirements of the 
positivity of DE entropy, the entropy bound and the GSL 
cannot be met simultaneously. In the authors used 
the future event horizon to study the phantom thermo- 
dynamics and found that the GSL could be respected if 
the phantom entropy is negative. The problem with the 
future event horizon is that for the universe with DE with 
equation of state w —1, the thermodynamical descrip- 
tion breaks down on the event horizon Q. Furthermore, 
the definitions of the event horizon temperature and en- 
tropy could be less certain than a guess. Even if we use 
the similar temperature and entropy definitions of the 
apparent horizon for the future event horizon, the first 
law of the thermodynamics was not satisfied 

In the above discussion we have concentrated ourselves 
on the DE with constant equation of state. To study the 
thermodynamics of a dynamic DE, we will use the gen- 
eralized Chaplygin gas (GCG) [l3| as an example. When 
the universe is dominated by the GGG, the entropies of 
the apparent horizon and the radiation read [loj 

Sa = Sao ^c^^ C^) 

5, = 5.o(^) (8) 

where = [-w^o + (l + Wco)ia/ao)~^^^+"'>]^/^^+'''> . The 
entropy for the GCG can be obtained by solving Eq. H]), 
which can be expressed as 



Sc Tco 

SrO TrO 



4(n-|-3) 
Scl + 



(H-u;eo)^/(i+")r!;o'^V+3, 



RD, 

a ^ ao. 



(9) 



where Sci and Sc2 are integration constants. To have 
Sc > 0, the parameter n must satisfy the condition 
—3 < n < 3(1 -I- a). Numerical results show that this 
condition is not enough. For example, if we choose 



Wco ~ —0.88 and a — 1.57, which are the best fitting 
values from observations 17] , we find that Sw is negative 
after MD when n — 2. At late times, a — > oo, Sc — *■ Sc2- 
For positive entropy, Sc will be greater than Sr at late 
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times since Sr — > 0. The range of n to keep Sc positive 
can be more confined by numerical calculation. Choos- 
ing appropriate n to ensure Sc to be positive, we have 
shown the numerical results in Fig. [5] on the evolution 
of entropies of GCG, radiation and the apparent hori- 
zon. When n < 0, Tc increases with the expansion of 
the universe and the numerical results show that Sc can 
be less than Sr during RD and MD eras if Tco/Tro is 
large enough. If n > 0, then Tc decreases as the universe 
expands and Sc increases faster than Sr during RD and 
MD eras. When n = 1, the GCG and the radiation tem- 
peratures evolve in the same way and Sc can be larger 
than Sr during MD as shown in Fig. [2l It is clear from 
Fig. [2]that compared to the apparent horizon entropy, Sr 
and Sc are negligible, thus the entropy bound can be pro- 
tected for the GCG case. In addition, the GSL can also 
be saved in the GCG case, since the total entropy evolves 
basically in the same way as the entropy of the apparent 
horizon. Though in the GCG dominated period, Sr de- 
creases as the universe expands, owing to its negligible 
value compared to the apparent horizon entropy, its de- 
crease can be overcome by the increase of the apparent 
horizon. 
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FIG. 2: The evolutions of Sc, Sr and Sa with Wco = -0.88 
and Q = 1.57. The dotted lines are for (Sc/Sao) x (Tco/Tro), 
from top to down are n = —2, n = —1, n — 1 and n — 2 
(Note that Sc is negative after MD when n — 2). The dash 
line is for St/Sao, and the solid line is for Sa/Sao- 

In summary, in this work we have investigated the ther- 
modynamical properties of the DE. In calculating the DE 
entropy we have considered the volume of the universe 
enveloped by the apparent horizon and assumed that the 
physical volume and the temperature are related. The 
apparent horizon is a good boundary for studying cos- 



mology, since on the apparent horizon there is the well- 
known correspondence between the first law of thermo- 
dynamics and the Einstein equation [l6l |. Furthermore, 
it has been found that the apparent horizon is good in 
keeping thermodynamical laws Q. Assuming that the 
temperature of the DE has the form T ^ a^", we have 
derived the evolution of the DE entropy. For the DE 
with constant equation of state w > — 1, we have found 
the appropriate range of n for keeping DE entropy to 
be positive, which is the requirement of the statistical 
understanding of the concept of entropy. In this range 
of n, the entropy bound and the GSL can also be pro- 
tected. The negative point is that the allowed range of n 
for giving physically acceptable DE entropy leads the DE 
temperature to increase as the universe expands, which 
is different from the behavior of the thermal tempera- 
ture that decreases along the expansion of the universe. 
This conflict could be overlooked since the DE temper- 
ature and the thermal temperature may have different 
dispersion relations [1, [13] , and it is not necessary that 
these two different temperatures behave accordingly. In 
the era of phantom domination, the GSL requires that 
the phantom entropy increases as the universe expands 
and the phantom temperature T satisfies the condition 
|1 + 3w\Th/2 > T > 0. Since the horizon entropy de- 
creases to zero as the universe expands, the holographic 
entropy bound will be violated if the phantom entropy 
is positive. For the phantom with constant equation of 
state u; < — 1, we found that there is no common range of 
n so that the positivity of the entropy, the entropy bound 
and the GSL can all be satisfied. The physical require- 
ment on the DE entropy does not favor the phantom with 
constant equation of state. We have also extended our in- 
vestigation to the dynamical DE by using the GCG as an 
example. We have found that by appropriately choosing 
parameters, we can have positive DE entropy, and mean- 
while we can protect the holographic entropy bound and 
the GSL. Within the allowed parameter range for phys- 
ically acceptable DE entropy, the DE temperature can 
decrease and it can even scale in the same way as the 
radiation temperature does as the universe expands. 
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